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Abstract 

Model dependence of the capture rates of the negative muon capture in deuterium is studied 
starting from potential models and the weak two-body meson exchange currents constructed in the 
tree approximation and also from an effective field theory. The tree one-boson exchange currents 
are derived from the hard pion chiral Lagrangians of the NAnptjai system. If constructed in con- 
junction with the one-boson exchange potentials, the capture rates can be calculated consistently. 
On the other hand, the effective field theory currents, constructed within the heavy baryon chiral 
perturbation theory, contain a low energy constant that cannot be extracted from data at the 
one-particle level nor determined from the first principles. Comparative analysis of the results for 
the doublet transition rate allows us to extract the constant d^. 

PACS numbers: 12.39.Fe; 21.45.Bc; 23.40.-s 
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I. INTRODUCTION 



The study of the weak interaction in the deuterium at low energies is an important topic. 
Since the structure of the nucleon-nucleon (NN) force is at present well understood and the 
weak vector one- and two-nucleon interactions and the weak axial one-nucleon interaction 
are also well known, it can provide the most reliable information on the structure of the space 
component of the weak axial meson exchange currents (MECs) . This component contributes 
notably into the capture rates for the reaction of the negative muon capture in deuterium 

+ d — y + n + n, (1.1) 

from the hyperfine fid states. 

In reaction f ll.ip . the stopped muons are captured from the doublet or quadruplet hyperfine 
states with total angular momentum F=l/2 or 3/2, respectively. The corresponding capture 
rates, A1/2 and A3/2, were calculated in the past by many authors (for the references, see 
[1,2]). Since A3/2 ~ 10 s~^ and A1/2 ~ 400 s~^, only the doublet capture rate is of practical 
interest. 

Until recently, the capture rates were calculated with the deuteron and neutron-neutron wave 
functions derived from the variety of nuclear potentials of the first generation and with the 
MECs constructed in the tree approximation. The systematic investigation of the structure 
of the operator of the tree weak axial MECs based on the chiral Lagrangians was started in 
Ref . [sj . These MECs were applied in calculations of the capture rates and neutron- neutron 
spectra for reaction (11. 1 p in Refs. 0-0] • The obtained values of the doublet capture rate 
Ai/2 were in the interval 398 s~^-416 s~^, which was considered in a reasonable agreement 
with the measured value A1/2 = 409 it 40 s~^ [sj, in contrast to another measurement 
providing somehow larger A1/2 = 470 ± 29 s~^. Let us note that in these calculations 
the MECs are of the one-boson exchange form. If used in conjunction with the nuclear 
wave functions generated from one-boson exchange potentials (OBEPs), all the parameters 
of the calculations are fixed. The same space component of the weak axial MECs induces 
other very important weak processes in the two-nucleon system, such as the weak neutrino- 
and antineutrino-deuteron disintegration and the proton-proton fusion. In the terrestrial 
conditions only the reaction (II. ip has the real hope to be studied experimentally with an 



accuracy of ~ 1 % in the near future [10|. This is challenge also for theorists to improve 
the calculations of the capture rates. Let us note that the above described approach was 
highly successful in describing the electroweak processes in nuclei at low and intermediate 



energies [ll|, ll2||. We shall call it the Tree Ap proximation Approach (TAA). It is also called 
the Standard Nuclear Physics Approach [isj. 

A more fundamental approach to the study of electroweak processes in nuclei was possible 
after the appearance of effective field theory (EFT) [13, |15| . In this way one constructs 
also terms of higher order, including loops besides the leading order of a current (potential) 
operator. The expansion is done in a perturbative parameter q/A^ << 1, where g is a 
quantity characterizing the system (momentum, energy,...), which is small in comparison 
with the heavy scale A^ 



Recently, the reaction (II. ip was considered in two varieties of such an EFT. In Ref. 16 
Ai/2 was calculated within the framework of a pionless EFT. In this approach, an unknown 
low energy constant (LEC) Li^^ enters the weak axial MECs. It cannot be extracted from 
the one-body processes nor derived from fundamental considerations. In our opinion, the 
application of the pionless EFT to the reaction (II. ip is doubtful, because in this process 
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the contribution from the pion pole (the induced pseudoscalar) is important. Since the pion 



is absent, Chen et al. [16[ introduced the pion pole by hands which is a misconception. 
Moreover, the heavy scale turns out to be of the order of the pion mass m^r, A^^, ~ m-,^ ~ 
140 MeV, whereas the value of the momentum transfer q is of the order of the muon mass 
m^, q ~ ~ 100 MeV. The calculated prediction for A1/2 for reasonable values of 
Li,A -5-6 (see FIG. 3 and TABLE I of Ref.[i6j) is ki/2 ~ 370 - 380 s'^ MeV, which is 
by ~ 5 - 10 % less than the results obtained within the above discussed TAA. 



In Ref. [13[, the muon capture in deuterium was studied within a 'hybrid' approach: the 
current operator was constructed within the framework of a heavy baryon EFT and the 
nuclear wave functions were obtained from the Schroedinger equation solved with the second 
generation Argonne fig potential [l^. The weak axial MECs operator of the heavy baryon 
EFT contains the one-pion exchange part and also a (contact) short-range term, the strength 
of which is given by the LEG rf^, that should be fixed by data. In Ref. [l3| . is taken 
from the analysis of the triton /3-decay rate 18| to predict the Ai/2(^>S'o) for the transition 



d — 7- ^5*0. Let us note that the value of d^ is correlated with the cutoff A entering a Gaussian 
regulator (strong form factor). Since the Argonne potential is not a OBEP containing such a 
form factor, the value of A is not fixed by the NN scattering data. Rather, the extraction of 
d^ is made for three values A = 500 MeV, 600 MeV and 800 MeV. For these values of d^ and 
A, the doublet capture rate calculated in Ref. [13] for the channel d — )• ^^o, Ai/2(^>S'o) = 245 
s~^, is by ~ 5.5 % less than the value Ai/2(^5'o) = 259 s~\ obtained in the TAA calculations 
0. 

Here we calculate the capture rates for the reaction (11. ip within both the TAA and the 
hybrid calculations. In both methods the basic symmetries of quantum chromo dynamics, 
such as the gauge symmetry and the spontaneously broken chiral symmetry, are reflected. 
The nuclear wave functions are obtained by solving the Schroedinger equation 

Hm = Em, H = T + V, (1.2) 

where T is the the kinetic energy operator and V is the NN potential. Since our TAA weak 
currents, constructed from chiral invariant Lagrangians in Refs. 0, 0, 0, |l^ 20|, are of the 



QBE type, we use OBEPs in Eq. (11.21) . too. Such potentials are successfully applied in the 
TAA for describing nuclear phenomena at low energies. This modeling of the nuclear force is 
a surrogate of calculations based directly on the quark-gluon dynamics because of the non- 
perturbative feature of the quantum chromodynamics in this region of energies. The well 
known first generation Bonn OBEPs [21I, 2^ contain the exchanges of scalar, pseudoscalar 



and vector bosons. In the BNN vertices, phenomenological form factors of the form 

^"(f^) - (lif )"" (1-3) 

are applied, with q the three-momentum transfer, Aq the so called cutoff mass and nQ=l 
or 2 depending on the specific couplings. These form factors reflect the extended structure 
of the nucleon. Later on, high quality second generation GD-Bonn OBEP appeared (23| . 
containing the BNN form factors of the type (II. 3p . too. Another type of OBEPs has been 
constructed by the Nijmegen group jl^]. These potentials contain the Gaussian strong form 
factors, 

FBNNif) = e-'^'/'^B , (1.4) 
where the cutoffs A^ are extracted from the fit to the NN scattering data. 
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In order to make our calculations of the capture rates consistent, we employ in Eq. (11. 2p 
the Nijmegen OBEPs and use the Gaussian form factors (11.41) and couplings entering these 
potentials both in our weak and EFT MECs. As it has been shown in Section 3 of Ref. [25| 
the potential current of the range B, j^^^(2,p.c.), of the weak vector nuclear MECs satisfies 
the nuclear Conserved Vector Current (CVC) equation 

g^j;^(2,p.c.) = [VB,m)], (1.5) 

where Vb is the OBEP of the same range B and Jq (l) one-body vector charge density. 

In its turn, an analogous current of the weak axial nuclear MECs, j^^ ^(2,p.c.), fulfills (see 
Section 3.1 and Appendix A of Ref. ^]) the nuclear Partially Conserved Axial Current 
(PCAC) equation of the form 

g,J5VB(2,p.c.) = [Vb, M^)] + ^^mlAUq')M%{2) , (1.6) 

where J5o(l) is the one-body axial charge density and A^^(2) is the associated two-nucleon 
pion absorption/production amplitude. Further, f^^ is the pion decay constant and is 
the pion propagator. So the potential current in this case depends only on the employed 
potential model. Our definition of the MEC operators guarantees that there is no double 
counting if the MEC effects are calculated with the wave functions from Eq. (II. 2p . 
Let us note that the presented approach has already been successfully applied to the descrip- 



tion of nuclear phenomena. For example, in Refs. [26|, |27[, the cross section of the reaction 
of the backward deuteron disintegration, 

e + d^e' + n + p, (1.7) 



is well described up to energies of 1 GeV j28|, l29| . In the weak sector, the capture rate of 
the reaction 

/i- + ^He ^ u^ + ^H, (1.8) 



was calculated in Ref. 30|. The obtained value Pn = 1502±32 s ^ is in a very good agreement 
with the result of the precise experiment HH ij, Pq'^p = 1494.0 ± 4.0 s'K 



The procedure of implementing phenomenological strong form factors into the electromag- 
netic MECs of the OBE type was in detail discussed in conjunction with the Bethe-Salpeter 
equation in Ref. [33|. Similar procedure for the weak axial MECs was performed in Ref. (33 |. 
The construction of our weak axial nuclear MECs in conjunction with the Schroedinger equa- 
tion (II. 2p stems in Ref. [20] from these currents. The set of Eqs. (II. 5p guarantees the gauge 
invariance of the calculations with the weak vector current, whereas the set of Eqs. (II. 6p 
provides the correct relation between the weak axial current and the related pion produc- 
tion/absorption amplitude. When one deals with potentials of particular form such as the 
OBEPs containing phenomenological form factors, the splitting of the continuity equation 
for the currents into the set of equations for particular exchanges follows naturally. On the 
other hand, in the given case and also for a potential V of general form one can generate 
conserved currents in the weak vector sector using the minimal substitution [35, For the 



weak axial sector, an effective procedure for constructing the MECs satisfying the nuclear 
form of the PCAC, 

q,r,^{2,p.c.) = [V, J- (1)] + ^^mlAUq')M^{2) , (1.9) 

is not known. 
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Since the most important weak axial MECs are the A(1232) excitation currents of the tt 
and p ranges, we pay particular attention to them. At first we generalize them so that also 
pieces depending explicitly on the off-shell parameters Y and Z |37| are included. Besides, 
we construct the A excitation currents also from the 'gauge symmetric' Lagrangians 38|-|40 . 
In calculations of the capture rates, we take specific values of the parameters Y and Z, which 
allows us to study the model dependence. In particular, if we take Y and Z as obtained 
from the requirement that the vriVA and pNA interactions should not change the number 
of degrees of freedom of the free A liil-iaj , we obtain the same capture rates as calculated 
with the A excitation currents constructed from the gauge symmetric Lagrangians. 



In our hybrid calculations, we apply the weak MECs that were used in [13|. However, we 
require that the weak vector MECs satisfy Eq. (11.51) and that the leading order terms of the 
weak axial MECs satisfy Eq. (11. 6p . also for the Gaussian strong form factor. In order that the 
weak axial MECs satisfy the nuclear form of the PCAC (11.61) . one needs a potential current 
constructed in [3]. In the weak vector sector, in addition to the pion pair and pion- in- flight 



terms, we employ the A excitation current of the vr and p ranges. These vector MECs are 



well known from the study of the process n + p^d + ^aX the threshold in the TAA j45 



Since the A excitation currents saturate about 30 % of the MEC effect, we find it necessary 
to include them here, too. Let us note that the strength of the A excitation current of 
the TT range constructed within the TAA from the gauge symmetric vriVA Lagrangian is 



close to the strength of such a current constructed within the heavy baryon EFT [46|. The 



inspection of Table 2 [46| also confirms that this current should contribute non-negligibly. 
As we shall see, our calculations confirm that it is really the case. 

In numerical calculations of the doublet capture rate, we use in the hybrid approach the 
same value of the cutoff A^ as in the TAA, which allows us to extract from the comparison 
of the doublet capture rates a unique value of the LEC d^. 

In Section [Til we discuss briefly the methods and inputs necessary for the calculations, in 
Section [TTll we display the equations for the capture rates and in Section HVt we present the 
results. We conclude in Section |Vl Further, in Appendix \^ we present the weak currents 
of the TAA and then we collect the Fourier transforms of the used MECs. In Appendix [B] 
we consider analogously the EFT currents, and in Appendix O we deliver all the multipoles 
of the currents used in the numerical calculations of the capture rates. 



II. METHODS AND INPUTS 

To obtain the capture rates one must calculate first of all the matrix elements of the weak 
nuclear currents between the initial and final nuclear states. Here we describe the needed 
ingredients of these calculations. 



A. Weak nuclear currents of the TAA 

The weak hadron current, triggering the reaction (11.11) . is 

= + h^- (2-1) 

The weak vector and weak axial nuclear currents j~ and j^, respectively, consist of the 
one- and two-nucleon parts presented in Appendix |A1 There is practically no uncertainty 
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associated with the one-body part. Hence we concentrate on the effects of the two-body 
currents in the channel d — )■ ^Sq. 

The weak axial nuclear MEC J5^(2) that we consider here is of the OBE-type with the tt 
and p exchanges. It can be divided [44| into the potential and non-potential currents. The 
potential current of the range B, j^_^(2,p.c.), satisfies the nuclear PCAC equation (11. 6p . 
The main part of the non-potential weak axial exchange currents contain the model in- 
dependent p-TT current and the A excitation currents that are model dependent. In our 
calculations, we shall adopt the vrA^A and pNA Lagrangians used for many years [s^, 47 
to study the vrA^ reactions and the pion photo- and electroproduction on a nucleon (model 
I) and also the gauge symmetric Lagrangians proposed recently [s^, 11] . In model I, we de- 
rived the A excitation MECs from Lagrangians possessing the hidden local SU{2)l'x. SU{2)ii 



symmetry [40[. In particular, the vertices containing the A isobar field were chosen as 



+ C 



■NAp ; 



where 



NAna 



nNA 



^>;fO^,{Z)m ■ {d,n + 2Ugpa,) + h. c. 



f 



NAp 



NAp 



G 

-gp^^^fO^^iX)^^^,'^ ■ prju + h. c. 
gp^^^>,fO,,{X){d,m)-pr,. + h.c.. 



(2.2) 

(2.3) 
(2.4) 
(2.5) 



Here T is the operator of the isospin 1/2 — )■ 3/2 transition. The constant fnNA is extracted 
from the A isobar width 48 with the result 



7. 



NA 



/47r = 0.767 ± 0.024 fm^ 



(2.6) 



Besides, Gi = 2.525 [37[. Further, the operator 0^y{B) is taken in the form |37l. 1471. |49|. |50 



G{B) 



6. 



p,u 



C{B)^,^, 



2 + ^ 



(2.7) 
(2.8) 



The parameters X, Y and Z do not influence the on-shell properties of the A isobar, hence 
they are called off-shell parameters. They were systematically extracted from the data in 



Refs. m. 47, 49, 50 



On the other hand, the parameters entering the Lagrangians (12.21) were restricted in 4ll, |42 



on the basis of field theoretical arguments by the requirement that they should conserve 
the same number of degrees of freedom of the A isobar as possessed by the free one. This 
resulted in that it should be 



y = 0, 



z = l 

2 



6-2 = 0. 



(2.9) 



This was criticized in [37[, because the requirement G2 = fixes the ratio of the multipole E2 
to multipole Ml amplitudes kinematically, which is unacceptable and the restriction (12. 9p 
was refused as a whole. In our opinion the consistency requirement G2 = means that the 
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model I can be only applied if the amplitudes derived from the Lagrangian C\^p^ Eq. (12.51) 
do not contribute. 

Later on, the problems related to the Lagrangians (12.21) were reconsidered in [3 



39, 43 



In particular it was shown in j43| within the framework of the Hamiltonian formalism with 
constraints that the result 2' = | persists. 

Since the contribution to the effects of the A excitation currents for the process (II. ip from 
the Lagrangian (12. 5p is negligible, we consider in our calculations two sets of parameters Y 
and Z. In the set la we take 

Y = Z = -\. (2.10) 

In its turn, the set lb is defined as 



y = 



(2.11) 



The A excitation currents corresponding to the set (I2.10p were used in all previous calcula- 
tions of the MECs effects for the reaction (II. ip . 



39| do not contain any off-shell pa- 



The vriVA and piVA Lagrangians of the model II [38|, 
rameters. They are constructed in such a way that they conserve the degrees of freedom of 
the free A isobar. We write them as 40 



'-NAp 



rgs. 

NAnai 



fn 



NA 



e^iuap [{dfi^u) Tlhla ^] ■ {dfiT^ + 2f^gpais) + h.c. 



(2.12) 



MMa 



(2.13) 

The values of the coupling constants are obtained from the condition that the new La- 
grangians, Eq. (I2.12P and Eq. (I2.13P and the standard Lagrangians, Eq. (12. 3 p and Eq. (12. 4p . 
respectively, are equivalent for the on-shell A isobar. 

It turns out that the A excitation currents of this model differ from those of the model I, 
set la (IXTOj) . only by the factor (M/Ma)^ ^ 0.58 [M(Ma) is the nucleon (A isobar) mass]. 



B. Weak nuclear currents of the EFT approach 

For the estimation of the MECs effect in the hybrid calculations, we use the MECs operator 
of Ref. [ill from which we omit the second term at the right hand side of Eq. (19). This 
term is suppressed by the factor 1/M in comparison with the leading term, which provides 
by itself only a small contribution to the weak axial MECs effect. In the weak vector sector, 
Ando et al. [l3[ present in Eq. (17) the non-relativistic version of the currents Jemi^) 
Jem{b) of the set (C.4), constructed within the heavy baryon EFT 46[. They coincide with 



our TT-pair term (1A9P and pion-in-flight term (lAlOp . constructed within the formalism of 
hard pion Lagrangians. Here we add to them the A excitation currents of the vr and p 
ranges (lAlip and (lA13p . It can be shown that the A excitation current of the vr range 
(lAlip . if taken with the set (I2.10p and multiplied by the factor (M/Ma)^ corresponds to 
the A excitation current Eq. (C.4) of Ref. j46|]. The TAA calculations show that 

the currents flAlip and flA13P contribute non-negligibly. So in order to assess the capture 
rate Ai/2{d — ?■ ^Sq) correctly, one should include explicitly also these currents in the EFT 
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calculation of the capture rates. In our opinion the current f lAlSP can be obtained from the 
heavy-fermion Lagrangians of Appendix C 46| as well. 

In the weak axial sector, we add to the weak axial MECs [isj the vr potential term f lA17p . 
This ensures that also the weak axial MECs of the EFT approach satisfy in the leading 
order the PCAC constraint ([US]) [ij. 

The weak axial MECs of the heavy baryon EFT contain the known LECs Cj, {i = 1, 2, 3, 4), 
Ce and the unknown LECs di and d2- These two constants enter in the calculations of the 



observables effectively in the combination di + 2d2 or as used in the formalism of Ref. [13 
in the combination 

= di + 2d2 + -C3 + ^C4 + ^ . (2.14) 

3 3d 



On the other hand the combination di + 2^2 is also expressed in terms of the LEC cd a.s [51[ 

di + 2d2 = -^^CD , (2.15) 
where MAr=0.939 GeV is the nucleon mass and A^ ^ 1 GeV. 

Let us note that the constant d^ (cd) enters not only the short range part of the weak axial 
MEC but also the pion production/absorption amplitude in the NN collisions and the NNN 
force. 

The dependence of the LECs on the off-shell parameter Z was discussed by Bernard et al. 
[5^ in the model where the LECs are saturated by the heavy mesons and the A resonance. 
In parallel with the model la fl2.10p . in which Z = — | we obtain the set Ila 

C2 = 3.37, C3 = -4.70, C4 = 3.31, (Ila), (2.16) 

whereas for the case lb 02.111) in which Z = }^ the LECs of the model lib are 

C2 = 1.98, C3 = -3.31, £4 = 2.61, (lib). (2.17) 

Another set of these constants, which we design as the set lie, has been extracted in js 
from the data, 

C2 = 1.67 ±0.09, C3 = -3.66 ±0.08, C4 = 2.11 ± 0.08 , (lie). (2.18) 

If one omits from the LECs C2, C3 and C4 the contribution from the A isobar, one gets the 
model Ild 

£2 = 0.047, C3 = -1.371, C4 = 1.63, (lid). (2.19) 



For the constants Ci and cg we take 13|, [53 



ci = -0.60 ±0.13, C6 = 3.70. (2.20) 
New set of the LECs has recently been delivered in Ref. [s^]. We take it as model He, 

ci = -0.85 , C2 = 3.1±0.2, C3 = -4.4+Jq, £4 = 3.3 , (He). (2.21) 

These new LECs are close to the set Ila f l2.16p obtained for Z = However, these new 
LECs are extracted with much larger errors than the set lie (12.181) . 
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C. Nuclear potentials 



We use the Nijmegen I (NI) and Nijmegen 93 (N93) [23] OBEPs. The couphngs and cutoffs, 
entering these potentials, are employed also in the MECs. In particular, the cutoffs are 

= 827.5 (NI) , A^ = 1177.11 (N93). (2.22) 

Let us note that our main results are related to the high quality second generation NI 
potential with x^/Nrfata=l-03, whereas for the N93 potential x^/Nrfata=l-87. 



III. CAPTURE RATES 



The capture rates A1/2 and A3/2 are related to the statistical capture rate, kstat, as 

1 2 

Ktat = 2 Ai/2 + -A3/2. (3.1) 



Using the method of Ref. [55|, we write the statistical capture rate for the process fll.ip in 
terms of the multipoles. 



A 



stat 



~3 



Gf cos 6c(f)f,{0) 



+ \<Xjf\\Lj-Mj\\d>f}, 



(3.2) 



where 



A 



M„(A - - z/V4M„) 



2 Mr, 



1 + A/M„) 



(3.3) 



Mn {md) is the neutron (deuteron) mass, and = —0.00267 MeV is the binding energy of 



the muon. Further in ([32]), the weak interaction constant Gf = 1.16637 x 10"^ GeV"^ |48| 
cos 6'c=0. 9749 and the wave function of the bound muon at the center of the deuteron is 
0/x(O) = (^/i,r tt)^''^/v^, m^^r is the reduced mass of the muon and a is the fine structure 
constant. 



The multipoles and their reduced matrix elements are defined in Section 4.1 of Ref. [20 



Generally, the capture rates Kp for the hyperfine state F can be obtained from the equation 

Kf = Ktat + Cf5K, (3.4) 

where 

J i F^ 



J,; 1 



(3.5) 



and (5A depends on the nuclear dynamics. Here Jj is the total angular momentum of the 
initial nucleus. 

For the reaction (II. ip . Ji=l and 



1 

3' 



1 

6' 



(3.6) 
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whereas 

SA = V6Mn 



Gp cos 9cM^) 



TT 

J J' 1 







Xjf,JJ' 



J' J.. 1 



^^^[-1 1 



1 1 Qj<\jf\m'-Tr'\\d><\jf\\tTp-TjT'\\d>* 

' < XJf\\^ff ~fr^\\d>< Xjf\\Lj, ~ Mj,\\d>' 



+ ( o' J ) <>^3f\\Lj-Mj\\d><\jf\\Lj,-Mj,\\d>* 



(3.7) 



Here J = \/2J + 1 and the symbols ( "'^ "'^ ) and { "'.^ "'.^ "'.^ i are Wigner's 3jm and 



nil rn2 rn^ 



6j symbols, respectively |56l |. 



J4 J5 J6 



IV. RESULTS AND DISCUSSION 



Here we present first the results for the capture rates obtained for the reaction f ll.ip in the 
formalism of the TAA and then in the hybrid calculations. 



A. Results for the TAA 



In the formalism of the TAA, we calculated the contributions to the capture rates from all 
channels d — )■ '^''^'^^Lj^, where L=S,P,D,F, j j=0,l,2 and from the multipoles J=0, 1,2,3 of the 
one-nucleon currents. The contribution of the weak MECs was taken into account in the 
multipole J=l and in the channel d — )■ ^5*0. We also estimated the MEC effect due to the 
weak vector MECs j~{p.t.), Eq. (|A9|) . and j "(tttt), Eq. (lAlOp in the channels d -)■ ^-Po,i,2- 
The results of the calculations of the doublet capture rate A1/2 for the reaction fll.ip are 
presented in Table [H 

It is seen from Table [I] that the results depend on the potential and current model used. As 
noted above, our basic model is Nl/Ib: the potential NI qualitatively supersedes the N93 one 
and the current model lb should be preferred because the parameters of the A excitation 
currents are restricted by the reasonable demand that the number of degrees of freedom 
of the A isobar should be conserved. Moreover, the same results as displayed in the third 
row of Table [T] are obtained with the potential NI and with the current model based on the 
Lagrangians (12. 3p and (12. 4p . As is seen from the fifth column of the Tabled the resulting 
MECs effect for this case is 16.7 s~^, which is ~ 4 %. The estimated error in A1/2 for this 

model due to the 3 % variation of the constant (I^^ma^ Mtt given in Eq. (12. 6 p is ~ 1 s~^, 
which is ^ 0.25 %. 

Let us note that the time component of the weak axial MECs contributes by +2.3 (+2.9) s~^ 
for the models Nl/Ia and Nl/Ib (N93/Ia and N93/Ib). Our calculations show that the effect 
of the weak vector A excitation currents is ~ 4 - 6 s~^. It is a non-negligible contribution 
to Ai/2 of the order of 1 - 1.5 %, which cannot be neglected if the expected error of the data 
is ~ 1.5 %. If one looks for the effect of the heavy meson exchanges, one obtains for the 



10 



TABLE I: Partial contributions to A1/2 (in s"^). The value of the constant ('^^j /47r = 0.783 
fm^ is used. In the first column, the potential and the current model are displayed. In the second 
column, the contribution from the one-body current in the channel ci — ?• ^5*0 is given. In the third 
column, the contribution of the MECs in the channel d — t- ^Sq is presented, whereas in the fourth 
column, the contribution of the one-body current to all considered channels but d — t- ^Sq is given. 
In the fifth column, we give the total MECs effect and in the last column, all the contributions are 
summed up. 





lAo 


AMECo 


AIA MEC lA+MEC 


Nl/Ia 


239.2 


22.0 


160.4 23.8 


423.4 


Nl/Ib 


239.2 


14.8 


160.4 16.7 


416.3 


N93/Ia 


238.8 


29.1 


160.4 30.8 


430.0 


N93/Ib 


238.8 


22.0 


160.4 23.8 


423.0 



model Nl/Ib that 

AAf/f^ = 24.3 s-^ - 7.6 s-^ = 16.7 s'S (4.1) 

where the first number at the right-hand side of the equation is due to the MECs of the 
pion range and the second number stems from the heavy meson exchanges. So the effect of 
the short range MECs is ~ 30 % of the long range one. 

In Table [TTl we present the contribution from particular multipoles to the capture rates 
calculated with the lA currents. It is seen that the contributions from higher multipoles not 
taken into account cannot change the results much. 

TABLE II: Partial contributions to the capture rates (in s~^) from the multipoles J, calculated in 
the lA approximation. The neutron-neutron '^^'^^Lj^ partial waves with L=0,l,2,3 and jj=0,l,2 
are taken into account. In the second column, only the contribution from the channel d — t-^ Sq 
to the multipole J=l is considered. In the third column, the capture rates are calculated with 
the multipoles J=0 added; in the fourth column, all the multipoles with J=l but arising from the 
channel d — t-^ 5*0 are added; in the fifth (sixth) column, the multipoles J=2 (J=3) are added. 



J 


1 CSo) 12 3 


Astat 
Al/2 
A3/2 


83.2 93.2 119.6 139.4 141.0 
239.2 249.3 322.4 394.6 399.6 
5.1 15.2 18.3 11.9 11.7 



We now give for the model Nl/Ib the final results for the capture rates. In the channel 

Al, = 88.1s-\ A?/2 = 254.0 s-\ A^/, = 5.2 s-' . (4.2) 
The full calculations provide 

Astat = 146.4 s-S Al/2 = 416.3 s-\ A3/2 = 11.4 s'^ (4.3) 

11 



The result A1/2 = 416.3 s~^, Eq. (14. 3p . seems to be in agreement with the one of Eq. (38) of 
Ref. [gI], Ai/2 = 416±7s~^. However, it is more correct to compare this result with the value 
Ai/2 = 430 s~^ of our Tabled obtained from the first generation realistic potential N93 and 
the TAA current model la, because this sort of potentials and of the current model was used 
also in j^. In Ref. the value A1/2 ~ 400 s~^ was reported, as the result of calculations 
with similar potentials and currents as in [6j . On the other hand, according to Ref. 13|| , the 
contribution of the higher partial waves was later reevaluated 0] and an enhancement of the 
Ai/2 by ~ 10 s~^ was achieved. Then A1/2 ~ 410 s~^ is in good agreement with the result of 
[61], but it is by 5 % smaller than our corresponding value of A1/2 = 430s~^. Let us note the 
calculations of Ref. [s^] reporting A1/2 = 402 s~^, also performed with the first generation 
realistic potentials and the current model la. 



B. Results for the EFT currents 



Here we provide the results of calculations of the LECs and c/j by comparing the dou- 
blet transition rate A'j'/2 calculated with the weak MECs of Section IIV A[ We made the 
calculations with the potential NI and the weak MECs discussed in Section IIIBI Using 
the Goldberger-Treiman relation we connected the coupling of MECs [l3| with the constant 
gnNN that we took from the potential. Since the cutoff A^ is also taken from the poten- 



tial, our hybrid calculations are consistent as much as possible. As in Ref. [13[, we use here 
the weak form factors in the linear approximation in the expansion in the four momentum 
transfer g^, given in Eqs. (IBip . and (]B2|) . However, this provides only a small difference in 
the results, in comparison with the full dependence. 

It follows from Table [T] for the model Nl/Ia that A'j'/2 = 261.2 s^^. This provides for the 

model Nl/IIa the value of the LEC = 3.225. In both models, Z = For the other 
considered models, we take = 254.0 s~^ obtained for the model Nl/Ib and the resulting 

constants d^ and cd are presented in Table ITnl 



TABLE III: Values of the LECs and cd obtained by comparing the doublet capture rate for the 
channel d — )• ^Sq, calculated in Section [IV Al for the model Nl/Ib, and the one calculated with the 
weak MECs of Section [II B[ The value Z = 1 is taken in the model Nl/Ib and also in the model 
Nl/IIb, in calculating the contribution of the A resonance to the LECs Cj, i = 2, 3, 4. In the model 
Nl/IId, the contribution of the A resonance to the LECs Cj, i = 2,3,4 is omitted, whereas the A 
excitation current of the pion range is explicitly taken into account. 





Nl/IIb NI/IIc Nl/IId Nl/IIe 


d^ 

CD 


2.410 2.155 0.625 2.680 
2.173 2.436 -0.231 2.407 



Comparing the second and the fourth columns of Table IIIII shows that taking into account 
the A resonance effect by the method of the resonance saturation of the LECs and taking 
it into account explicitly by calculating the MECs effect is not equivalent. Also comparing 
the third and the last columns one finds about 20 % change in the value of the constant 
d^. Having in mind the large uncertainty in the LECs of the set He (I2.2ip one concludes 
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that the change is not essentiaL Let us also note that the value d^=2.68 was obtained with 
y = Z = — I in the A excitation currents. If one employs Y = and Z = the value of 

(i^=2.66, so it changes insignificantly. The value of the entering the NI potential f l2.22p 
is close to the value of one of the cutoffs, A=800 MeV, used in the analysis of the triton 
beta decay [l8[. For this value of the cutoff, the extracted = 3.90 ± 0.10 [ll], which is 
enhanced at least by 30 % in comparison with from our Table lllli 

We have also calculated the influence of the weak axial potential current j?;^^{2,p.c.) on 
the value of the constant d^. As discussed in Ref. [4^, this current is usually absent in 
calculations of the weak processes. If we omit it the value of the constant (i^=2.410 in the 
second column of Table [Till increases to ^"^=2.740, thus it changes by ~ 13 %. It follows 
that if one would like to extract the value of d^ with an accuracy better than 10 %, then 
one should take the contribution of the potential current into account. The calculations also 
show that omitting the potential current causes an enhancement of the doublet transition 
rate A1/2 by ^ 1 %. 

The constant c^ has recently been extracted [5l|, together with another LEG c^;, entering 
the contact part of the NNN force, from the data on the triton beta decay, binding energies 
and point-proton radii of the 3N system and ^He nucleus, with resulting value cd=-0.2. As 
it is seen from our Table UTTt only the value cd=-0.231, corresponding to the model Nl/IId, 



is in agreement with the analysis of Ref. j5l|]- Let us note that with the choice C4 = —63=3. 4 



5l| one obtains the value ^"^=1.1 using C£)=-0.2. 



It follows from our calculations that the effect of the time component of the weak axial 
MECs is ~ -1 s~^, which is in agreement with Ref. [3]. This is in contrast to the TAA 
calculations based on the hard pion Lagrangians, where the time component contributes as 
f» +2 s~^. The short range part of the hard pion time component reduces the value d^= 
2.155 (see the third column of Table UTTj) to d^= 1.91, if one uses this component instead of 
the soft pion one in fitting d^. 

Our = 254.0 s~^ differs by 9 s"^ from the same quantity, rjj° = 245 s~\ given in Table 
1 of Ref. [l3|. The main part of the difference can be assigned to the difference of ~ 7 s~^ in 
the contribution from the one-body currents, obtained from comparing IAo=239.2 s~^ (see 



our Table H]) with the value of r^rf=232 s ^ of Table 2 131. O n the other hand, our MECs 



effects, AMECo=14.8 s , are close to 13 s obtained in [13 



V. CONCLUSIONS 

We have evaluated the capture rates for the reaction of muon capture in deuterium (11. ip . 
both using the TAA currents and those derived within the EFT approach. 
The weak TAA currents, presented in Appendix [XJ are of the one-boson exchange type 
obtained from the hard pion chiral Lagrangians and they satisfy the nuclear CVC and 
PCAC constraints, Eqs. (11. 5p and (II. 6p . respectively. The final state neutron- neutron wave 
functions were generated from the high quality second generation potential NI and from 
the realistic potential N93 Since the potentials are also of the one-boson exchange 

type, employing in the TAA currents the same couplings and strong form factors (II. 4p . we 
performed fully consistent calculations, presented in Section IIV A[ For the main object of 
interest, the doublet capture rate for the channel d — J-^ Sq, we predict (see Table [H 
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model Nl/Ib) 

A°/2 = 254 ± 3s"^ (5.1) 

This result was obtained with the neutron-neutron wave functions derived from the NI 
potential. The error reflects the uncertainty in the vriVA and pNA couplings, possible 
effects of the neglected short range effects and applied approximations. In the model Nl/Ib, 
the ttNA and pNA couplings preserve the physical degrees of freedom of the free A isobar. 
The lA currents contribute to the value (15.11) of Ay 2 by 239 s~^, whereas the MEC effect is 
15 s~^, which is ~ 6 %. In the full calculations, we considered the neutron-neutron "^^^^Lj^ 
partial waves with L=0, 1,2,3 and j/=0,l,2, and the contributions of the lA currents to the 
multipoles J=0,l,2,3. For the full doublet capture rate we got 

Ai/2 =416 ± 6s-^ (5.2) 

In addition, the estimated error includes also the uncertainty due to the neglect of the 
contribution from the higher multipoles. 

The EFT currents that we used are discussed in Appendix [Bl The hybrid calculations of 
the capture rates accomplished with these currents are presented in Section IIVBI These 
calculations are consistent to the extent that we again use in the MECs the couplings and 
strong form factors from the potential NI. We extract the unknown LEG by comparing 
the capture rate A'j'y2 with its numerical value calculated with the TAA currents. Besides, 
we adopt various available sets of the known LECs Cj, discussed in Section [11 B[ As is seen 
from Table UTTl the value of changes within 25 % for various sets of Cj. The exception is 
provided by the case in which the contribution of the A is eliminated from q, and the A 
excitation current is taken into account explicitly. Then the value of d^ is suppressed by 
the factor 4. 

Comparing our results with those of Ref. [l3| we see that our calculations provide the value 



of A*!* 12 = 254 s ^ which is by ~ 4 % larger than the analogous value T^^ ° = 245 s ^ of 13 



by ^7 



'0. 



10 of the reac- 



Equally, our total capture rate A1/2 = 416 s~ differs from F^^ = 386 s~ [13 
In conclusion we stress that the planned precise experimental investigation 
tion (II. ip is of fundamental importance. It will stimulate efforts to understand better the 
details and limits of application of both the TA and EFT approaches and will certainly shed 
more light on the value of the important LEG d^ (cd)- 
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Appendix A: The weak currents of the TAA 



The hadron currents consists of the one- and two-nucleon parts. The one-nucleon currents 
are of the form, 



3 



9p 



a 



8M2 



[P "^a - {a- P)P + (a ■ q)q- i{P x q)] 



9AFAiq 

2M 



a-q)q 



[a-P) 



T 



9p 



2Mm, 



(Al) 
(A2) 

(A3) 
(A4) 



Here P = p' + p, q^ = p'^^ — p n, where p'^ {p^) is the four- momentum of the nucleon in the 
final (initial) state and the induced pseudoscalar form factor is 



gp{q^) = 2MgAm^AUq') ■ 

For the other weak form factors we employ the dipole parametrization, 

F^{q^) = 1/(1 + gVM^) , = 0.711 GeV' , 

FA{q') = 1/(1 + q'/Ml), Ml = 1.04 GeV^ 



We use for the constant gA the value [48 

9a = - 



1.2694 ± 0.0028. 



(A5) 



(A6) 
(A7) 



(A8) 



1. The weak exchange currents 

The two-nucleon part also consists of the weak vector and weak axial vector parts. We 
present first the weak vector MECs. They are 



1. The TT-pair term, 



fnNN \ ' 



J \p.t.) = F^, iq') A^(gt) i^.V(^2 ) ^1 (^2 ■ gl) ^ {n x r,r + (1^2). 

(A9) 

2. The pion- in- flight term. 



ttNN 



Fi{q^)qi {Si ■ qi) (0^2 ■ 92) 



qi - ^2 



AUq2)F^NA<l2 



t in X r2T + (1^2). 



(AlO) 
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3. The A excitation current of the vr range 



J/ (A) 



—I 



+ 



9 (Ma - M) 



1 + f{Y,Z) 



Q2 ^2 



1 - 2f{Y,Z) 



xAUq2') F^NN^ (^2 ■ gl) + (1^2) 



(All) 



where q = g7|g| and 



f{Y,Z) = (1-M/Ma) [C(r)+C(Z) + 2C(r)C(Z)(2 + M/MA)] 



C(a) = -i^ + a), Cr = 2Gi 



J'ttNA f; 



ttNN 



Mml 



(A12) 



4. The A excitation current of the p range 
i'^(A) = 



^^^^ F,V)gx^4 



9 (Ma - M) 



1 - 2/(r,F) 



(ga X (ct2 X gs)) 



+ 



1 + 4/(r,F) 



i{ai X (g2 X (a2 x ga))) « {n x fs)" 



xA^(g-'22)F^^^(g-^) + (1^2) 



where 



2M V M / ■ 



(A13) 



(A14) 



The current j '^(vrvr) of Eq. ( lAlOl) is written in such a form 25|, |33| that the potential current 

r{p-c.) =rip.t.) + /'^(vrvr) (A15) 

satisfies the CVC equation (11 .Sp for any form factor -Fttatat, if the one-pion exchange potential 
also contains it. Here we use for the Gaussian form factor (11.41) the following approximation 
in the current (lAlOp . 



^^^A-(g-;^)A-(g-^)F^V(?2) 

HI ~ H2 



oo / 



-1)" 



q2) + {<i2 +0 

n~l ■ 



n=l 



-•9 ->2 

q{ - q2 



Al 



^Fiq/)F^NNiq2 



1 - 



qi - qj 

2KI 



(A16) 



The weak axial MECs are 



1. The vr potential current [44 

' fTTNN\'^ 9 A 



2M 



FA{q^) [ (g + X Pi) + (Pi + lai x q) i (fi x fs)" ; 



xA^(g-'2^)P,V(C)(^2-g2) + (1^2) 



(A17) 
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2. The p-TT current, 



mP^) = -{^) ^[l+"^pA^(9t)][i'i + (l + <)^(^ixgi)] 



xF.N^iq,^) Aliq^^) F^NN{q2 ) (c?2 ■ gl) ^ x fa)" + (1 O 2) 

xFp^7v(gi') A^(g2') F^^7v(g2') (^2 ■ gl) ^ (ri x + (1 o 2) .(A18) 
Only the second part of Eq. (lAlSP contributes to the rate A 1/2 sensibly. 



3. The A excitation current of the pion range, 

9 a f f: 



2 



9(Ma - M) V 

x[l - gA-(g2)g-'.]|4 

+ [1 + f{Z, Z)] I (<Ti X gl) z (ri x f^T } 

X A^(gt) i^.';viv(g2') ('Js ■ gl) + (1^2). (A19) 



1 - \f{Z,Z) 



4. The A excitation current of the p meson range, 

X [1 - gA^(g2) g. ] { 4 [1 + /(F, Z)] x (0^2 x ^2) 

+ [1 - 2/(r, Z)] zai X (g2 x (0^2 x gl)) ^ (n x fa)" } 

X A^(g2') Fl^^{qi) + (1 O 2) , (A20) 

where 

5. The potential current of the p range 13 

h>-c-) = i^i)' (^Mf 9AFA{q') {r2" [g x {a, x gl) + x (A x {a, x gl)) 

+i(fi X [A X ((T2 X g2) + iai x (g x (0^2 x gs))]} 
xA^(g-'2^)F^V(g-'2^) 

- (I) ' 72^ ^ ^"^ • ^"^^ + ^(^^ " ^^^'^ ^^^^ ■ " 

X A^(g-'22) ^^^^(gl,^) + (1^2), (A22) 

6. The time component of the weak axial MEG 

P5W = -(^)'^[l+"^;A^(g-2)]F,^;v(g-l^)A^(g-^)F.^;v(g'2^) 

X (a2 ■ g2) i (fi X f2)'^ + (1 o 2) . (A23) 

In Eqs. flA9P - flA23p . Pj = p/ + Pi, g = gl + g2, qi = vl — Vi- In reaction (11. ip . the isospin 
components of the currents and charge densities O ~ = O ^ — iO are effective. 
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2. The Fourier transform of the weak MECs 



Here we provide the Fourier transform of the weak MECs in the same order as they are 
hsted in the previous section. We start with presenting the form factors arising due to the 
form factors of the type fll.4p after the Fourier transformation of the weak MECs containing 
one boson propagator. They are 
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(A24) 



(A25) 



/ AbXb rriB\ 
^"^y 2mB Ab) 



(A26) 
(A27) 
(A28) 



In Eqs. ( ]A24p -( |A26p . the function erfc{x) is the complementary error function 58 
First follow the weak vector MECs: 



1. The vr-pair term, 

fip.t.) = -^iJfIiF^{q')a,{a,■f)^{nxT2re^^'^■'^^W,^ + (1 o 2) . (A29) 
47r 

2. The pion- in- flight term, 



1 fU 



27r3/2 , 



q \ 



xE^^/llW + (1^2) 

i=l 



(A30) 



where 



O3 = (^1 • g)(a2 ■ V,-)V,-, 04 = 2(<Ji ■ V,-)(a2 • V,-)V, 



(A31) 



and 



L 



(A32) 
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M \ 2AI 



ho^ ■ (A33) 



Here /3 = (p^ + + m^)/2pq and Ql{P) is the Legendre polynomial of the second sort 



1581. 



Numerically, only the current (tttt) contributes non-negligibly to the A1/2 for the 
transition c/ — > IS'o, 



1 E 



TT 



L N=L±1 

K I J ]\h g I 



X [e ® [[ai ® aj^ ® [V^g) ® rA'(r)]^]^]^° + (1 o 2) , (A34) 
where e = e±, cq are the orthogonal unit vectors, and 

cl+1 = yrrr, cl_i = yZ, ^jv+i = VA^ + l, ^tv-i = VVV. (ASS) 

3. The A excitation current of the vr range, 

h (A) = -^ ge^Ml-M) ^^"^^^'^" E, ^""-^ ^^"^^^ + (1^2), (A36) 

d^(A) = a2r,\ (A37) 

d2"(A) = -f(a2-r)r2", (A38) 

dl{A) = 2(^1 X (T2) ^ (n X 72)'^ , (A39) 

OKA) = z(f X ai)(a2 • r) z {n x Ts)'^ , (A40) 

Fj; = — IVi. [1 - 2/(r, Z)] , = W2^ [1 - 2/(r, Z)] . (A41) 

4. The A excitation current of the p range. 



Jp(A) = -^ , J,; Fr(g^)e-fa--^) E x ^"(A) + (1^2) ,(A42) 



X p 

-V _ /TT/ 2 ^ 



= (w^2p - - w^ip) [1 + 4/(r, r)] , F/^ = w^2p [1 + 4/(r, r)] . (A43) 

Next we present the Fourier transform of the weak axial MECs: 
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1. The TT potential current, 

f2 

J : 



Otip.c. 
OUp-c. 

dt{p-c. 

F7n{p.C. 



2. The p-TT current. 



■kNN ^tt 

47r 2M 



10 



QAFAiq^) e^^^-^^^^ Yl F^AP■^■) ^Up-c) +(1^2), (A44) 



■x ai) (a2 ■ r) , d2"^(p.c.) = 1i (<T2 ■ r)(^i x Vi) r^" , (A45) 



%q{o2 ■ r) % (fi X T2T , Ol^{p.c.) = 2((72 ■ f) Vi ? (n x fs)" 
^g(c?2 ■ r) , d6",(p.c.) = (ai x q){a2 ■ f) % (ri x f^f , 

X CT2) r2" , 08^(p.c.) = -02 % (fi X T2Y , 
i((Ji X f)((T2 ■ r) , d;'o^(p.c.) = f{a2 ■ r) i (n x ra)" , 
F2^{p.c.) = F^^{p.c.) = Fi^{p.c.) = F^^{p.c.) = Fg^{p.c.) 

Fs-^{p.c.) = Win/xT, , 

FlOn{P-C.) = W2n ■ 



ttNN 



1 1 



3 v\ 



J AMqa 7rg V 27r -v ■ • p 

xz(ti X r2)"[Ai|5"(p7r) + As^Iptt)] + (1^2) 



Here 



L ^ N=L,L+2 jk 



X < 



1 L + 1 
1 k L 
1 j 1 



[C(i, fc, Ar)P° + E ^iv 



X /^IatW E (-i)'i^ < 

jk 



1 N L-1 
1 A; L 
1 i 1 



N=L-2,L 



[C{j.KN)\ 



where 



A2 

Trp 



1 1 

+ 



A2 ' A2 A2 A2 ' 

p TTp TT p 



(A46) 
(A47) 
(A48) 
(A49) 

(A50) 
(A51) 
(A52) 



(A53) 



(A54) 



(A55) 



ul = VL + 1, aL+2 - VL + 2, 6l_2 = - 1 , 6l = \/Z,(A56) 

(A57) 



FL{p,q) 



Qiia) - Ql{/3) 



p2 + a(l — a) + (1 — a) + a w? ' 



a = (p^ + + m^) / (2apg) , 
/[2(a - . 



[p^ + (1 - a)^g^ + m 



(A58) 
(A59) 



22 



Further the symbols [C{j, k, A^)]^'' are defined as 

k, iV)]^o = [e ® [Y^if) ® [a, ® a^]^] Y° , (A60) 

and the piece A2j^ (pn) can be obtained from the term Aijg (pvr), Eq. ( lA54p . by the 
change 

ml^l/a, FL{p,q) ^ QL{a). (A61) 

ai bi Ci 

At last, the symbols ^ 02 &2 C2 \ in Eq. (lA54p are Wigner's 9j symbols [5 

03 h C3 

3. The A excitation current of the ir range, 

•^^"^^ 367r(MA -M)\ J 

X [1 - qAUq'm E (A) + (1^2), (A62) 

^{t = F^t = 4W^2.[l-^/(^,^)], 

Fgt = -W^Al + fiZ, Z)] , Ft = [1 + fiZ, Z)] . (A63) 

4. The A excitation current of the p range, 

^'"^^^ ~ 367r(MA - M) ^^^^ 

X [1 - qA^q^m jz Ff^ dt{A) + (1^2), (A64) 

i=l 

Ff, = AiW2,--W,,)[l + fiY,Z)],F,^^ = AW2Al + fiY,Z)], 

Xp 

K = {W2p--W,,)[l-2f{Y,Z)], Ftp = W2p[l-2f{Y,Z)]. (A65) 

X p 

5. The potential current of the p range 

xOlpip.c.) + (1 2), (A66) 
dlp{p.c.) = g(<?i X f)(a2 • q) r,^ , O^iv-c) = -2z(ai x f)(a2 ■ Vi) r^'^ , (A67) 
Olip.c.) = X f)(a2 ■ r) r^'^ , Ol^ip.c.) = -iqf^a^ ■ q) i x r^f , (A68) 

6lp{p.c) = -2f(a2 ■ Vi) t (fi X , 6lp{p.c.) = -f{a, ■ f) t (n x f^f , (A69) 
Fip{p.c.) = F2pip.c.)=F4p{p.c.) = F,p{p.c.) = Wip, (A70) 
Fspip.c.) = F^pip.c.) = mpW2p . (A71) 
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6. The time component of the weak axial MEC 



~a, ^ ( I^NnV 1 a -^{l-a)+ia{q-f)+i{q-r^) 

X z (n X r^riA^ptipTr) + A^ptipn)] + (1 o 2) . (A72) 

Here 

Aip^M = ^ t-'^LYLoir) HI^{t) . (A73) 

L 

The function H'^^ir) is defined in Eq. (]A57p and the piece A2P5(pvr) can be obtained 
from the term Aip5"(p7r), Eq. (1A73I) . by the change ( 1A61I) . 

Appendix B: The EFT currents 

For the one-body currents, we take the currents from Eos. flAlD - (lA4l) of Appendix |Al with 
the form factors in the quadratic radius approximation [59|, l60[ , 

Fiiq^) ~ 1 - r?. = 0.59 fm% (Bl) 

6 

FA{q^) ^ 1 - ir^g^ r\ = (0.403 ± 0.030) fm^ (B2) 
6 

However, this approximation changes the results only a little bit in comparison with the 
dipole form factors used in the TAA calculations. 

1. The weak exchange currents 



For the weak vector currents, we take the currents from Eqs. flA9l) - flA13P of Appendix [XI So 



we add to the vr-pair and pion-in-flight terms considered in [13] the A excitation currents of 



the TT and p ranges. Inspecting Table 2 of Ref. |46| one can in addition expect a non-negligible 



contribution also from the pcuvr and two-pion exchange currents. 

As to the weak axial MEC operator, we adopt here the main part of this current used in 



13| and add to it the vr potential current of Eq. ( 1A17I) d emanded by the PCAC constraint 



(II. 6p . In our notation, the currents (19)-(21) of Ref. [13| are: 

^Ib = -^A^(gt) ■ q2) ^ in X nT + (1^2), (B3) 

9aFa 



A. 



2B 



2M/2 



1 - 1 

- Pi i (fi X fa)" + 2c3 q2 - (64 + -)i(ai x q^) i {n x fa)" 



+ ^i^i y< q) iinx fa)" 



^f{Q2) {^2 ■ q2) 

+ di{ai + aa t^) - ^2 «(^i x (Xa) i (fi x fa)"} + (1 o 2) , (B4) 



^(^"i' = -T^JVTf r,^ AW2 ) (^2 ■ q2) + (1^2). (B5) 



2M/2 

In Eq. flB3l) . we keep only the leading order of the time component considered in Eq. (18) [l3| . 
because the contribution of this part of Aa^ is small and the correction to it is suppressed by 
the factor f» 1/M^. Let us note that the time component flB3P is the soft pion approximation 



to its hard pion form (IA23p . 
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2. The Fourier transform of the EFT weak axial MECs 

Multiplying the currents (IB3|1 . flB4|) and (IB5|) by the form factor squared of the type (11 ■4p 
we obtain 



Xb = -^^^-^^^{^2 ■ r)e^^'-'''^ W,^ I (ri x f^f + (1^2), 
47r 4/^ 



(B6) 



A. 



2B 



i^^^'^^'"^ ^ ^^^^^^ ^^"^^^^ +(1^2), (B7) 

dl{2B) = iq{a2 ■ r) % x ra)" , 02^(25) = 2 (0^2 • r) Vi 2 (n x fs)'' , (B8) 

03^(25) = -a2 2 (ri X , 04^(25) = f(a2 ■ r) z (fi x r^T , (B9) 

d«(2S) = 2a2T,\ dl{2B) = -2f{a2-r)T,\ (BIO) 

0^(25) = -i{ai X ^2) ^ (ri X fa)" , (5^(25) = z(<Ji x f)(a2 ■ r) i (fi x fa)" , (Bll) 

d<'^{2B) = -(ai X g-)(c?2 ■ f) ^ (fi X f2)" , d?o(25) = <Ti r^" , (B12) 

Fi(25) = F2{2B) = W^jAm^, F^{2B) = WijAx^, F^{2B) = W2JA, 

F,{2B) = c^W.Jx^ + d,Wj2, F^{2B) = c^W^. , 

Fj{2B) = {c, + ^)W,Jx^ + d2W^, Fs{2B) = (64 + ^)1^2., 

Fg{2B) = ^^^Wijm^, Fio{2B) = d,W^ . (B13) 

hqu 92) = ^^^^^ (a2 ■ r)e^^^""^^ W,^ r,^ + (1^2). (B14) 
47r M/^ 

Appendix C: The multipoles of the currents 

We first present the contribution to the multipoles from the lA currents. In order to make 
the equations more transparent, we do not write the argument gr/2 in the Bessel functions 
ji{qr/2), unless the argument differs, which is the case of the p-vr current. The factor a/2 
arising from the isovector matrix elements is not kept in the reduced matrix elements of 
the current, but is included in the overall constants in front of the integrals in Eqs. ( 13. 2p 
and (13. 7p . We also take into account the factor l/^n, entering the reduced matrix elements 
according to the definition given in Eq. (4.20) of Ref. |2o|, by keeping l/ztg in the integration 
volume in the same Eqs. (13. 2 p and (13. 7p . 

1. The multipoles of the lA currents 

a. J=0 multipoles 
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+u,{r)/V2) +iV2^-j/ (I I " drul,^,iK,r) {(jo + J2) K(r) 
-Uo{r)/r) + [(jo + J2)«2W + (2jo - 5j2/2)M2(r)/r] V2} 



+00 



+ / dru\,^(K,r)ji -u'o(r)/r + uo(r)/r'' + (u'2(r)/2 



+U2{r)/r) /V2r\ 



(CI) 



<=^Pl||M50|M> = 



2M m,, JO 



+00 



+00 



M 7o 
+2ii2 (r)/r) /V2 



dr o(/«, r) jo [u'{r) - Uo{r)/r + (M2(r) 



(C2) 



b. J=l multipoles 



<^ S,\\Tr'\\d >= dTul,^,{K,T) \joUo{T)-32U2{r)/V2] , (C3) 



V f+00 



-j2U2{r)/V2 



8M2- 



+00 



+ '^[7^ Jo ^'<oA-.r) 



iiKW - + ^J^ dru'^Qo^^{K,r) [{7ji 

-3js) u'^{r) + (14ji + 9j3)M2(r)/r] 

X f+00 

+UQ{r)/r^)\ --J^ dr 4o,i{'^, r) [{jo - j2/2)u'^{r) 
+Sjou'2{r)/r + 3j2U2{r)y] } , (C4) 



So\\L,,\\d> = -z[gAFA{l 



-»2 

gpq 



2Mm, 



q 



f+00 
■/o 



2M2(,?^ 



[(ji + j3)M2(r) + (2ji - 3j3)u2{r)/r] 

2 r+00 

- o / ^00,1(1^^ ^) [Oo + j2X(r) + 3^2 (-'Uo(r)/r 

\i \/2 /■+°° 
+uo{r)/r^)\ - dru^oo^^{K, r) [(jo + j2)M2 W 

+3jo'U2(r)/r - 6j2?/2(r)/r^] } , 



(C5) 
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<' So\\M5i\\d> = 



. ? / p 9pQo^ 

It^KQaFa , 

2M JQ 



+00 



■ QaFa 



+00 



+2xi2(r)/r)] , 



(C6) 



<3por57"'ii(i> 

<3Po||Mi|M> 



<=^Pl||f57"^||(i> 



M '2 



3 



M rl 



V 2 2 

, Pr\\Mi\\d>=V2>F^Jl, 



(C7) 
(C8) 



(C9) 



<=^Pi||Li||d> 
P2\\iTf\\d> 

<^P2\\Ml\\d> 



<^ P'2|KT5llM> 
<^ P'2||M5l||c?> 



M 
V 2M 



IL+'^F^K,, 



, <' P2\\frn\d>^ ^9 aFaIL, (CIO) 



qG 



V 

M t1 



V5. 



-VbF^'jl,, ^P,\\L,\\d>^j^Fr 



-2P,^i^,\ 



<^P2||Li||d> = W--^ 



3 P^ [3^,1 ^1 



5 M 



8A 



2M 



M 7-1 
'5 > 



• 9 / ^ofi'Pwi • /7;9aFaji 



<^ D2\\L^i\\d> = -tV2{gAFA 



Q 9p wi 
2Mm,, ^ ' 



(Cll) 
(C12) 



(C13) 
(C14) 

(C15) 
(C16) 

(C17) 



II = drM°i_2(«,^)jiN(r) + M2(r)/\/2], 



^2^ 



4 



(C18) 
(C19) 

(C20) 
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r+oo 



d^ul,^{K,r)ji U2{r) 



r+oo 

/ di- uloA'^^ ^) b^uoir) - (2jo + j2)M2(r)/V2] , 

^0 

r+oo 

/ drulQ^{K,r) [j2Uo{r) + (jo - J2)u2ir)/ V2] , 
Jo 

/ drulo JK,r)ji Dluo{r) - V2Dlu2{r) , 
Jo ' ^ J 



/I 



7I 



r+oo 

\ dru^i o{'^,r)ji [uo{r) - V2u2{r)] , 
Jo 

/ dru\,2iK,r)ji [uo{r) +U2{r)/V2] , 
Jo 

r+oo \ 

druli^^{K,r)ji [uo{r) - -^M^)] , 
drul,^{K,r)jiU2{r) , 



r+oo 



+9j2D 

+00 





-9j2Dl 



Jo 

K2(r)} , 

dru\^2i.i^,r) 

r+co 

/ dru^,2{K,r) 
Jo 

+932DI] U2{t)} , 



+00 







dru[,2{n,r) 



-9j2D', 

f+00 



U2[r 



•)}. 

r+oo 

druli^^{K,r) 
+9j2Dl]u2{r)} , 

r+oo 

/ dru^i^{K,r) 
Jo 

-9j2Dl\ U2{r)} , 

r+oo 

druli^^{K,r) 

r+oo 

/ dru^i^{K,r) 
Jo 



{3o + 32)DUo{r) 



5V2 



V2 



(10jo+j2)^i 



(jo - 2j2)Dluo{r) - Y [(5jo - J2)D\ 



{jo + j2)Dluo{r) + 



IOV2 



(10jo + j2)I^i 



(jo - 2j2)i^^«o(r) + ^ [(5jo - j2)Dl 



{jo + j2)^>o(r) - ^ [(lOjo + j2)Di 



{jo - 2j2)D\uo{r) 



25^2 



(5jo-j2)^i 



joi^i - 2j2 



dr 



U2{r) , 
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J=2 multipoles 



<' D2\\iff\\d> 



2M 



mag i 
52 I 



(C38) 



^2 



<^ Pi\\L52\\d > = i{gAFA- 

.QaFa 



^ )ll<'P,\\M,,\\d>- 



2Mm„ 



2M 



{qaFa 



3 

72^ 



2Mm„ 



<3p2||T'^'||rf>-^ 



2M 



<^P2||M52|M> 
<'F2\\lfi\\d> 



^ r2 I T?v 7-2 



^V3^(^.i^.-^)/?A + ^V3 



^ r2 



9 G^I^x 



+6V2F^'ll 



13A 



<3F2||L52||ci> = iV2{gAFA 



q 9p 

2Mm„ 



)/|,,<3F2||M52|M> = 



/2 , ,-,/2^^/2 



(C40) 



(C41) 
(C42) 

(C43) 
(C44) 



/■ + 0O 

/ c?'^^20,i(«^'0i2K('^) +'fi2(r)/v2] , 

^0 

drM}i_2('«>^)bi^*o(r) + v2(-2ji + 3j3)M2(r)/5] , 



'3A 
r2 



+00 



— 

-'4,A - 



r - 



-'7A — 



/■ + 0O 

/ drii^i,;,(«,r)[j3Xio(r) + (-9ji + ldjs)u2{r)/10V2] , 
Jo 

f dru\^2{i^:r)[ji'^o{r) - (4ji + 9js)u2{r)/5\/2] , 
Jo 

r+oo 

/ dr uI^^^{k, r) 32 [i?>o(r) + D\u2{r)/^/2\ , 

J 

/■ + 0O 3 



(C45) 
(C46) 
(C47) 
(C48) 
(C49) 
(C50) 
(C51) 
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-'SA — 


r+oo 

Jo 


"31, aI'^' 






-'9A — 


r+oo 

Jo 


"ll.Al"'' 


n,^(r'] -1- 71^ ii^(r\ /a/91 


\\jOO j 


-'lOA — 


r+oo 

Jo 




r\ \ iir^(r\ -I- ii^(r\ / \/9] 

1 ) J2 L-f^+tto V ) T u_a2\i J/ V , 




-'ii — 


r+oo 

Jo 




f)32 U2{r)/r , 


\\jOO ) 


-'12A 


r+oo 

Jo 




r)j2 U2{r)/r , 


(C56) 


-'13A — 


r+oo 

Jo 




r)j2 U2{r)/r. 


(C57) 



d. J=3 multipoles 



<3p2r57^iM> 



<'F2\\f^r\\d> 



<^D2\\fr'\\d> 
<^ D2\\Uz\\d> 



'lA 
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P2\\Ls\\d>-- 



^3 



pV 7-3 



35V5M 



2JlgAFAll,,<'F2\\tff\\d>=^ 



1 7-3 



F^llx 



-3 r II f lU ^_ pVj^ 

7M ' 



, <^ F2\\L,\\d>-- 



= -2i 



.qG 



V 

M 7-3 



2M 



7|, <ir'2||iT|i||d>=-2i^^FA7|, 



(CSS) 



(C59) 
(C60) 



%9p_ 



2Mm 

+ ^^^^^ ' 



2M 



(C61) 



-'lA 



-'2A 



In 



-'5 



-'6A 



druli^^{K,r)j3U2{r), 

p+oo 

/ c^'^Wsi A('^'^)i3K('^) - U2{r)/5V2] , 

JO 

(•+00 g 

/ drulo^{K, r)\j2Uo{r) - V2(j2 + 7:j4)u2{r)/7] , 
Jo ' 2 

/ dru^Q AK,r)\j2UQ{r) - V2(j2 - 6j4)M2(r)/7] , 
Jo 

P + OO 

drulo^,{K,r)UDluoir) - V2Dlu2{r)] , 

r+oo 

dr uI^^^{k, t) [UJ2DI + {j2 + 15u)Dl]u2{r) , 



(C62) 
(C63) 
(C64) 
(C65) 
(C66) 
(C67) 
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-'7A 



'8A 



-'9A 



f+oo 

f ^ 4 



(J2 + ^J4)/^^ 



+ 2(3j2 + ^j4)I^i 



/■+00 c 4 ^ 

n2(r 



4 
25 



(J2 - ^J4)i?^ 



+ ^(3j2-^j4)/^i 



(C68) 



(C69) 



(C70) 



We now present the multipoles J=l of the TAA MECs given in Appendices lA II and lA 2[ 



2. Multipoles J=l of the TAA MECs 

First follow the multipoles of the weak vector MECs. 
1. The TT-pair term j (p.t.), Eq. (jA29]): 

V2 



So\\Tr'\\d> 



IttNN ^1 , 
IT Jo 



+00 



dr Uqq I 

(C71) 



2. The pion-in- flight term j (vrvr), Eq. (|A30|) : 



<i s'o||rr^lM> 



r2 _ f2 ^ 
00 -'^20J 



+J2(i^l2-i^02) Mr) - j2iF^o-Fio) 



+ {jo + 2j2){F^2-F22) 



U2{r 



(C72) 



3. The A excitation current of the vr range (A), Eq. ( 1A36I) : 



^1 s^\\Tr'\\d> 



97r(MA-M) [Jo 



X j2Mo(0 - (jo + j2/2)y2M2(r) 

f+oo 

-f{Y,Z) / drM°oi('^,0(W^2.-3iyi^/xO 

JO 

x[-2joUo(r) + v^j2M2(r)]} , 



(C73) 



4. The A excitation current of the p range j (A), Eq. (IA42P : 



97r(MA-M) [Jo 
32Uo{r) - (jo + j2/2)v^U2(r) 
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+8f{Y,Y) rfrM°o,i(«:,r) {W^p - ?>Wip/xp) 
xboWo(r)-j2M2(r)/v^]} , (C74) 



Now we write down the multipoles of the weak axial MECs, 

:^ a 

1. The TT potential current i^^{p.c.), Eq. (1A44I) : 



-2j2Dl) W^Jm^+j2W2-,]u^{r) + {[gji/2 
-(2jo+j2)W^2j2}v^M2(r)) , (C75) 

+4j2Dl) Wijm^ - 2j2W2^]uo{r) + {[-gji 

+ {32-3o)W2^]V2u2{r)) , (C76) 



2. The p-TT current {pT!'), Eq. (IA53P : 



<i5o|Kr5l||rf> = e-'^^(^-")/^^p[2/7 + J7] , (C77) 

<'5o||^5ilM> 



6y27r2Mc/Ag 

TTi. J Gnrng^q ^ ' ' ^ ^ 



n = / dru',,^,{K, r)3,{hqr)[{Hl, - Hl,)u,{r) + (if2'2 - Hl^)u2{r) / ^2] , (C79) 

JO 

_ r+oo p , 

/7 = dr «°o,i(«:, r) J2{hqr) [{Hl^ - //^^(r) + {-2Hl, - 

+HlJb + IQHlJl + 108/72'4/35) M2(r)/v^] , (C80) 
h = 1/2 - a. (C81) 

These multipole contributions correspond to the part proportional to Aijg (pvr) of the 

current (lA53p . The contributions due to the part ^2]^ (p^r) can be obtained from 
Eqs. flCTTl) and (IC78|) by the change (IMB . 



3. The A excitation current of the vr range j5^(A), Eq. ( 1A62P : 
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X rfrM°o,i(K,r) W2n[j2Uoir) - (jo + J2/2) M2(r)] 

+-fiZ,Z) / drul,iK,r){W2^-3W^Jx^) 
2 Jo 

X [-2joMo(r) + V2j2U2{r)]} , (C82) 



1 f+00 

--fiZ,Z) / cirMOo,(/«,r)(iy2.-3iyiJx^) 

Z JO 
x[joMo(?^) + V2j2^i2(r)]} • 



(C83) 



4. The A excitation current of the p range j^^^A), Eq. (IA64P : 



<'So\mi\\d> 



97r(MA - M) 



+00 



X / dru^oQ^^{K,r) W2p[j2Uoir) - {jo + j2/2)V2u2{r)] 



So\\U,\\d> 



-AfiY,Z) / drul^,iK,r)iW2,-3W,p/xp) 

<j 

x[joMo(r) - j2M2(r)/V2]} , 



(C84) 



97r(MA - M) 



+00 



c^^^ioo,i('*'^) ^2p[-j2Mo(?^) + U2 - 3o)u2{r) / \/2] 



+00 



-2/(F,Z) / drul^^^{K,r){W2p-m^p/xp) 

•J 

X joMo('^) + V2j2U2{r) } . 
5. The potential current of the p range j^p{p.c.), Eq. (jA66|): 

r+00 p -, 



(C85) 



+ 





+00 



dr4o,i{^^,r) {[g(8ji + 3j3)/5 - 2{2jo + j2)Dl]Wi,/m, 



So\\L,,\\d> 



+ {2jo+j2)W2p}u2{r)/V2) . 
. /(7p\Vm,\3 (1+^ 



(C86) 



V\2 



+00 



drulQ .^{K, r) {qji + 2j2Dl)Wip/mp - j2W^2p Uo{r) 
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+ dru',,^,{K,r) {[g(3j3 - 7ji)/5 + 2{jo - j2)Di]W,p/m, 
+ {j2-Jo)W2p}u2ir)/V2) . (C87) 

6. The time component of the weak axial MEC p'^{pTx), Eq. ( ]A72p : 



+ [j2{bqr)H^,o + Uoibqr) - j2{bqr))H^,^]V2u2{r)} 

+ - drul,^,{K,r)3^{bqr) {« + 2Hl,)u,{r) 

-[Hi, + {Hi, - ^Hl^mV2u2{r)]}) . (C88) 

This muhipole contribution corresponds to the part proportional to Aip'^{p7T) of the 
current f lA72p . The contribution due to the part A2P5(p7r) can be obtained from 
Eq. fOSj) by the change flMTjl . 



Next follow the multipoles J=l of the EFT MECs presented in Appendices IB II and IB 2[ 
3. Multipoles J=l of the EFT MECs 

~ 

1. The time component of the weak axial MECs, A2B, Eq- (IB6p . 

So\m,\\d>=-t^^^^ dru',,^,iK,r)j,W,Auoir)-V2u2{r)]. (C89) 
4vr/^ Jo 

2. The A2B term, Eq. (iBTl) . 

<'So\mi\\d> = i^^^ dTul,^,{K,T) ({-(3/4 + 2c4 + C3)joW^i./a;. 
V27rM/2 Jo 

+ [(3/4 + 2C4 + C3)jo + (C3 - C4)j2]W^2./3 + [-2(jo + ]2)D\/?, 

+ (1 + C6)gji] Wi^/Am^ - {di + 2d2)joW^} uo{r) 

+ {(3/4 + 2C4 + C3)j2Wij2x^ + [(C4 - C3)jo - (3/8 + C4/2 

+C3) J2] W2J3 + [2joDl/3+j2iDl + 9Dl)/15 

+ (1 + C6)gji/2] W^i./4m^ + (Ji + 2d2)j2Wj2} 

xV2u2{r)) , (C90) 



<i5o||L5ilM> = ^^^^^ dru1,,iK,r) {{-(3/A + 2c, + Cs)joW,^/x^ 
2TrMj^ Jo 
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+ [(3/4 + 2c4 + Cs)jo + 2(C4 - C3)j2]W^2 J3 + [2(2j2 - Jo)Dl/3 
+qji] Wi^/Am^ - (rfi + 2d2)joW^] uoir) 

+ {-(3/4 + 2C4 + C3)32Wi^/x^ + [(C4 - C3)jo 

+ (3/4 + 2c3 + C4)j2] W2J3 + [2jo/^i/3 - 2j2(/^i + 9D^)/15 
-gji/2] iyi./4m^ - (rfi + 2d2)j2W^j2} V2M2(r)) . (C91) 



3. The contribution due to the P{qi,q2) term, Eq. flB14p 



So\\L,^\\d>=^'-^^c,AUq^) 



+00 



dru^QQ jK,r)jiWi^ [uo{r) 



V2u2{r)]. 
(C92) 
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